Can the X(3872) be a 1 ++ four-quark state? 
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We use QCD spectral sum rules to test the nature of the meson X(3872), assumed to be an exotic 
four-quark (ccqq) state with J PC — For definiteness, we work with the current proposed 

recently by Maiani et al [J], at leading order in a s , consider the contributions of higher dimension 
condensates and keep terms which are linear in the light quark mass m q . We find Mx — (3925± 127) 
MeV which is compatible, within the errors, with the experimental candidate X(3872), while the 
SU(3) breaking-terms lead to an unusual mass-splitting Mx b — Mx = —(61 ± 30) MeV. The mass- 
difference between the neutral states due to isospin violation of about (2.6 ~ 3.9) MeV is smaller 
than the value (8±3) MeV proposed in [J. For the 6-quark, we predict Mx b = (10144±106) MeV for 
the Xhibbqq), which is much below the BB* threshold in contrast to the BB* molecule prediction Q, 
and for the X^(bbss), a mass-splitting Mx s —Mx b = —(121 ±182) MeV. Our analysis also indicates 
that the mass-splitting between the ground state and the radial excitation of about (225 ~ 250) MeV 
is much smaller than in the case of ordinary mesons and is (within the errors) flavour-independent. 
We also extract the decay constants, analogous to f n , of such mesons, which are useful for further 
studies of their leptonic and hadronic decay widths. The uncertainties of our estimates are mainly 
due to the ones from the c and b quark masses. 

PACS numbers: 11.55.Hx, 12.38.Lg , 12.39.-x 



In August 2003, Belle reported evidence for a new narrow state in the decay B + — » X(3872)if + — > 
J/ip7r + n~ K + 3], which has been confirmed by three other experiments Jj]. The A(3872) is the best 
studied of the new cc-associated states, A(3872), A(3940), F(4260), etc. [J It has a mass of 3872 MeV 
and a very narrow width T < 2.3 MeV at 95%. Upon discovery, A" (3872) seemed a likely candidate for 
?/> 2 ( 3 -D 2 ) or tp3( 3 D 3 ) [5], but the expected radiative transitions to \c states have never been seen. The 
7T7T mass spectrum favors high dipion masses, suggesting a J '/ip p decay that is incompatible with the 
identification of A(3872) — > 7r + 7r~ JM as the strong decay of a pure isoscalar state. Belle's observa- 
tion of the decay A (3872) — > J/ipj p determines C = +, opposite to the charge- conjugation of the 
leading charmonium candidates. The same paper [6( also reports the observation of the X decaying to 
J/ip ir + TT~ir°, with a rate which is comparable to that of the J/-07r + 7r~ mode. This decay suggests an 
appreciable transition rate to J/ip u) and establishes sizeable isospin violating effects. Finally, an analysis 
of angular distributions supports the assignment J PC = 1 ++ , but the mass of A (3872) is too low to be 
gracefully identified with the 2 3 Pi charmonium state. More recently, the Belle collaboration reported a 
peak in D°D°ir° which can be interpreted as the dominant decay mode of the X [7|. 

The anomalous nature of the A has led to many speculations: tetraquark 01 Hi , cusp [§| , hybrid [Toj . 
or glueball [HI]. Another explanation is that the A(3872) is a DD* bound state [H, [u|TlI G3, El] , as 
predicted before its discovery. 

In this work we use QCD spectral sum rules (QSSR) (the Borel/Laplace Sum Rules (LSR) [13, EE El 
and Finite Energy Sum Rules (FESR) [H, [H, ) , to study the two-point functions of the axial vector 
meson, A(3872), assumed to be a four-quark state. In previous calculations, the Sum Rule (SR) approach 
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was used to study the light scalar mesons [22], |23j, |24j, |25| and the Dj"j(2317) meson [2g,|27|, considered as 
four-quark states and a good agreement with the experimental masses was obtained. However, the tests 
were not decisive as the usual quark-antiquark assignments also provide predictions consistent with data 
and more importantly with chiral symmetry expectations p^ . I23L [H, |2^|. In the four-quark scenario, 
scalar mesons can be considered as S-wave bound states of diquark-antidiquark pairs, where the diquark 
was taken to be a spin zero color anti-triplet. Here we follow ref. and consider the X (3872) as the 
J PC = 1 ++ state with the symmetric spin distribution: [cg]g=i [cg]s=o + [c<z]s=o[c<7]s=i- Therefore, the 
corresponding lowest-dimension interpolating operator for describing X q is given by: 

U = ieab ^ ec [( g JC7 5 c b )(g d 7 M Ccf ) + (ql C7^c 6 )(g d 7 5 CcJ)] , (1) 

where a, b, c, ... are color indices, C is the charge conjugation matrix and q denotes a u or d quark. 

In general, one should consider all possible combinations of different 1 ++ four-quark operators, similar 
to e.g. done in [3l[ for the ++ light mesons and consider their mixing under renormalizations [12] 
from which one can form renormalization group invariant (RGI) physical currents. However, we might 
expect that, working with a particular choice of current given above will provide a general feature of the 
four-quark model predictions for the A (3872), provided that we can work with quantities less affected by 
radiative corrections and where the OPE converges quite well 1 As pointed out in [![, isospin forbidden 
decays are possible if X is not a pure isospin state. Pure isospin states are: 

X(I = 0) = ^±^L, and X(/ = l) = ^-^. (2) 

If the physical states are just the mass eigenstates X u or Xd, maximal isospin violations are possible. 
Deviations from these two ideal situations are described by a mixing angle between X u and Xd 

Xi = X u cos 9 + XdSinO, 

X h = -X u sin 6 + X d cos 9. (3) 

In ref. [![, by considering the X decays into two and three pions, a mixing angle 9 ~ 20° is deduced 
and a mass difference 

m(X h )-m(Xi) = (8±3)MeV. (4) 

In this work, we want to test in which conditions the results of the sum rules are compatible with the 
above predictions. 



II. THE QCD EXPRESSION OF THE TWO-POINT CORRELATOR 

The SR are constructed from the two-point correlation function 

IV(g) = i Jd 4 x e^<0|T[^)jt(0)]|0) = -H^ 2 )^ - M^) + n (q 2 )^. (5) 

Since the axial vector current is not conserved, the two functions, Hi and Ho, appearing in Eq. (O are 
independent and have respectively the quantum numbers of the spin 1 and mesons. 

The fundamental assumption of the sum rules approach is the principle of duality. Specifically, we 
assume that there is an interval over which the correlation function may be equivalently described at both 
the quark and the hadron levels. Therefore, on the one hand, we calculate the correlation function at the 
quark level in terms of quark and gluon fields. On the other hand, the correlation function is calculated 
at the hadronic level introducing hadron characteristics such as masses and coupling constants. At the 



In the well-known case of baryon sum rules, a simplest choice of operator 33] and a more general choice _34] have been 
given in the literature. Though technically apparently different, mainly for the region of convergence of the OPE, the 
two choices of interpolating currents have provided the same predictions for the proton mass and mixed condensate but 
only differs for values of higher dimension four-quark condensates. 
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quark level, the complex structure of the QCD vacuum leads us to employ the Wilson's operator product 
expansion (OPE). The calculation of the phenomenological side proceeds by inserting intermediate states 
for the meson X. Parametrizing the coupling of the axial vector meson 1 ++ , X, to the current, in 
Eq. (fTJ) in terms of the meson decay constant fx as: 



{Q%\X) = V2f x M x e, 



the phenomenological side of Eq. © can be written as 

yrphenf 2\ _ ^fx^X 



Ml 



(6) 



(7) 



where the Lorentz structure projects out the 1 ++ state. The dots denote higher axial-vector resonance 
contributions that will be parametrized, as usual, through the introduction of a continuum threshold 
parameter sq. 

In the OPE side, we work at leading order in a s and consider the contributions of condensates up to 
dimension eight. We keep the term which is linear in the light-quark mass m q , in order to estimate the 
mass difference in Eq. Keeping the charm-quark mass finite, we use the momentum-space expression 
for the charm-quark propagator. The light-quark part of the correlation function is calculated in the 
coordinate-space, and then Fourier transformed to the momentum space in D dimensions. The resulting 
light-quark part is combined with the charm-quark part before it is dimensionally regularized at D = 4. 

The correlation function, IIi, in the OPE side can be written as a dispersion relation: 



n 



OPE 



ds 



p(s) 



(8) 



where the spectral density is given by the imaginary part of the correlation function: ^p{s) = 
lm[Uf PE (s)}. After making an inverse-Laplace (or Borel) transform of both sides, and transferring the 
continuum contribution to the OPE side, the sum rule for the axial vector meson X up to dimension-eight 
condensates can be written as 2 : 

fS 



2f x M x e 



-M-* x /M 2 _ 



s/M A 



int.':. 



P (s) +nr 



<<?<?> 



(M 



where 



with 



p(s) = p pert {s) + p m «(s) + p^(s) + p {G2) {s) + P mix {s) + p^\s) 



(9) 
(10) 



p pert {s) = 
p m *(s) = - 



p<™\s) 



2 3 7T 4 



da 



or 

da 
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d/3 



/3)(1 + a + p) [{a + (3)m 2 c - a(3s 



(qq) [m 2 c - a(l - a)s] 
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T 
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m 2 c (qq) 
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2 5 7T 4 
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2 5 7r 2 a/3 2 
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d/3 



-3 ^(1 



A, 

l-Q 



da 



d(3 



p 2 [(a + P)m 



a + P) [{a + P)m 2 c - aps] , 

m 2 (l- (a + P) 2 ) (1- 2a -2/3) 



aPs\ 



2a 



[(a + P)m 2 c - aPs] 



(11) 



2 We have not included the effects of a dimension 2 term induced by the UV renormalon, [35l which we expect to be 
numerically negligible like in the other channels [37] , though this result needs to be checked. Instanton-like contributions 
which appear as a high-dimension operators will also be neglected like some other higher dimension condensate effects. 
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where the integration limits are given by a m i n = (1 — yl — 4rnf/s)/2, a maa: = (1 + yl — 4to 2 /s)/2 and 
{ftmin = am 2 ) / (sa — m 2 ) . We have also included the dominant contributions from the dimension-five 
condensates: 



p mix (s) 



m c (qga.Gq) 



l-Q 



(ml - a(l - a)s) + j d{3 [(a + (3)m 2 c 

0n 



a(3s\ 



a + f3 



(12) 



where the contribution of dimension-six condensates (g 3 G 3 ) is neglected, since assumed to be suppressed 
by the loop factor 1/167T 2 . The usual estimate (g 3 G 3 ) ~ lGeV 2 (a s G 2 ) [13] would deserve to be checked 
in more detail. We have included the contribution of the dimension-six four-quark condensate: 



(s) 



m 2 c (qq) 2 
12tt 2 



4to? 



(13) 



and (for completeness) a part of the dimension-8 condensate contributions 



3. 



n7 



(M 2 ) = - 



m 2 c (qga.Gq)(qq) 
24tt 2 



da 



1 



a(l-a)M 2 2(1 -a) 



exp 



a(l - a)M 2 



(14) 



III. LSR PREDICTIONS OF M x 




FIG. 1: The OPE convergence in the region 1.6 < M 2 < 2.8 GeV 2 for s\ /2 = 4.17 GeV. We start with the 
perturbative contribution (plus a very small m q contribution) and each subsequent line represents the addition 
of one extra condensate dimension in the expansion. 



In order to extract the mass Mx without worrying about the value of the decay constant fx , we take 
the derivative of Eq. ([9]) with respect to 1/M 2 , divide the result by Eq. ((9]) and obtain: 

M , _ i;z^ 



We should note that a complete evaluation of these contributions require more involved analysis including a non-trivial 
choice of the factorization assumption basis |38j . We wish that we can perform this analysis in the future. 
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This quantity has the advantage to be less sensitive to the perturbative radiative corrections than the 
individual moments. Therefore, we expect that our results obtained to leading order in a s will be quite 
accurate. 

In the numerical analysis of the sum rules, the values used for the quark masses and condensates are 
(see e.g. Ql [M III 53) 4 : 



m c (m c ) = (1.23 ±0.05) GeV, 
m u = 2.3 MeV, 

Tfiq = (m u + m d )/2 = 4.3 MeV, 
m s = 100 MeV, 
(qga.Gq) = m 2 (qq) with 
( g 2 G 2 ) = 0i 



m b {m b ) = (4.24 ± 0.06) GeV, 
m d = 6.4 MeV, 
(qq) = -(0.23±0.03) 3 GeV 3 , 
(ss)/(qq) =0.8 ±0.2, 



0.8 GeV , 



GeV 4 . 



1/2 



4.1 GeV, si' 2 



(16) 



4.2 



We evaluate the sum rules in the range 2.0 < M 2 < 2.8 for two values of sq: 
GeV. 

Comparing the relative contribution of each term in Eqs. (fTTj) to (|14[) . to the right hand side of Eq. 
we obtain a quite good OPE convergence for M 2 > 1.9 GeV 2 , as can be seen in Fig. [TJ This analysis 
allows us to determine the lower limit constraint for M 2 in the sum rules window. This figure also 
shows that, although there is a change of sign between dimension-six and dimension-eight condensates 
contributions, the contribution of the latter being smaller, where, we have assumed, in Fig. [1] to Fig. 01 
the validity of the vacuum saturation for these condensates. The relatively small contribution of the 
dimension-eight condensates may justify the validity of our approximation, unlike in the case of the 
5-quark current correlator, as noticed in [42l ]. However, the partial compensation of these two terms 
indicate the sensitivity of the central value of the mass prediction on the way the OPE is truncated. 



4.2 r 




2.3 2.4 

M 2 (GeV 2 



FIG. 2: The OPE convergence for Mx in the region 1.9 < M < 2.8 GeV for s 



1/2 



4.1 GeV. We start 



with the perturbative contribution plus a very small m q contribution (long-dashed line) and each subsequent line 
represents the addition of one extra condensate dimension in the expansion: + {qq) (solid line), +(g' 2 G 2 } +m q (qq) 
(dotted-line in top of the solid line), +{qq} 2 (dashed line), -\-m 2 i {qq) (solid line with triangles). 



In Fig. [2] are shown the contributions of the individual condensates to Mx obtained from Eq. (fT5|) . 
From Fig. [21 it appears that the results oscillates around the perturbative result, and that the results 



4 To leading order approximation in a s , at which we are working, we do not consistently consider the running scale 
dependence of these parameters. We shall use here the values of the quark masses obtained within the same QCD 
spectral sum rules methods compiled in |19| . They are defined in the MS-scheme, and have obtained within the same 
truncation of the QCD series from different channels and by different authors. 



G 



obtained up to dimension-5 are very close to the ones obtained up to dimension-8. For dcfinitcncss, the 
value of Mx obtained by including the dimcnsion-5 mixed condensate will be considered as the final 
prediction from the LSR, and the effects of the higher condensates as the error due to the truncation of 
the OPE. 



s =4.2GeV 
— Continuum 




1 — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 

2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 

M 2 (GeV 2 ) 



FIG. 3: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution. 

We get an upper limit constraint for M 2 by imposing the rigorous constraint that the QCD continuum 
contribution should be smaller than the pole contribution 5 . The maximum value of M 2 for which this 
constraint is satisfied depends on the value of sq. The comparison between pole and continuum contri- 
butions for Sq^ 2 = 4.2 GeV is shown in Fig. [3] The same analysis for the other value of the continuum 
threshold gives M 2 < 2.2 GeV 2 for 4 /2 = 4.1 GeV. 

4.3-, 

s 1/2 = 4.1GeV 

4 - 2 " s 1/2 = 4.2GeV 

4.1- 




3.8- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 

M 2 (GeV 2 ) 



5 More restrictive conditions are sometimes imposed in the literature, where, for example, it is required that the continuum 
contribution is smaller than 30 % of the total contribution. In this case no sum rule- window is allowed. In our analysis, we 
use a less restrictive criterion, having in mind that the role of the continuum is expected to be larger for high-dimensional 
current operators than in the usual p-meson channel, as indicated by different sum rules analyses in the existing literature. 
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FIG. 4: The X meson mass as a function of the sum rule parameter (M 2 ) for different values of the continuum 
threshold: sj /2 = 4.1 GeV (solid line) and Sq /2 = 4.2 GeV (dashed line). The arrows indicate the region allowed 
for the sum rules: the lower limit (cut below 2.0 GeV 2 ) is given by OPE convergence requirement and the upper 
limit by the dominance of the QCD pole contribution. 

In Fig.lU we show the X meson mass obtained from Eq. (|15p . in the relevant sum rules window, with 
the upper and lower validity limits indicated. From Fig. [4] we see that the results are reasonably stable 
as a function of M 2 . In our numerical analysis, we shall then consider the range of M 2 values from 2.0 
GeV 2 until the one allowed by the sum rule window criteria as can be deduced from Fig. [1] for each value 
the so— range of values. 

Using the QCD parameters in Eq. (fl6|) . we obtain the LSR predictions for different values of s and 
including the dimension-5 condensates: 

M x = (3908 ± 26 ± 13 ± 100 ± 46) MeV for sj /2 = 4.1 GeV, 

(3943 ± 30 ± 10 ± 80 ± 48) MeV for s^ 2 = 4.2 GeV. (17) 

The errors are due respectively to M 2 , (qq), m c and the truncation of the OPE. We have estimated 
the absolute value of the last error by varying the dimension-six and eight condensates from their vacuum 
saturation values to the ones where a violation of the factorization assumption by a factor two is assumed. 
The errors due to other parameters are negligible. One can notice that the main error comes from the 
uncertainties in the determination of the charm quark mass, which plays a crucial role in the analysis 
like in the one of other heavy quark systems. One can also notice that the central value of the mass 
prediction increases with sq. Apart the intuitive observation from an extrapolation of the known mass 
splittings from ordinary mesons which may not be applied for the multi-quark states (see e.g.(43j), So 
remains a free parameter. We shall try to fix its value using FESR. 



IV. FESR PREDICTION FOR M x 

As an alternative, we use the FESR, which can be obtained from Eq. ^ by taking the limit 1/M 2 — > 
and equating the same power in 1/M 2 in the two sides of the sum rules. We get up to dimension-six 
condensates: 



Equating the coefficients of the polynomial in 1/M 2 in both sides of Eq. (|18p gives n equations: 

2f x M%M% 1 = f ° ds s n p(s), n = 0, 1, 2... (19) 

J Ami 

Finally, dividing two subsequent equations (with n and n+1), we can obtain the mass Mx for any chosen 
value of n (which, formally, is expected to by the same for any n): 

C\ds s n+1 p{s) 

In contrast to the previous method, the FESR have the advantage of giving correlations between the 
mass and the continuum threshold so, which can be used to avoid inconsistencies in the determination 
of these parameters. Ideally, one looks at a minimum in the function M x (so), which would provide a 
good criteria for fixing both so and Mx- The results for different values of n are very similar, therefore, 
in Fig. O we only show the result for n = 0. One can see in Fig. [5]that there is no stability in sq, which 
can indicate the important role of the QCD continuum in the analysis. 

One can also notice, from Fig. 03 that the FESR converges faster than the LSR due mainly to the fact 
that here we use an expansion in 1/sq where ~ 4.1 GeV, while in the LSR, the expansion is done in 
1/Af 2 , where M ~ 1.4 GeV is much smaller. 




FIG. 5: The FESR results in Eq. (|20p for Mx as a function of so for n = including condensates up to different 
dimensions. The LSR results have been also inserted for easy comparison. 



V. FINAL QSSR PREDICTIONS FOR M x AND fx 

In order to exploit the complementary role of the LSR and FESR, we also show in Fig. [5] , the result 
obtained from the Laplace sum rule using M 2 = 2.2 GeV 2 . The factor K was introduced to account for 
deviations of the factorization hypothesis for the D = 6 [U HI, 52, El EE El and 8 [U H EH, El] 
—condensates; K = 1(2) refers to an assumption that respects (violates by a factor 2) vacuum saturation. 

The intersection point fixes the range of values of Sq to be: 

sjj/ 2 = (4.15 ± 0.03) GeV , (21) 

which is smaller than intuitively expected. This small value of the continuum threshold relative to the 
value of the resonance mass signals again the important role of the QCD continuum in the analysis, which 
is expected for correlators described by high-dimension current operators. Using this range of values of 
so, one can, definitely, fix the X mass to be: 

M x = (3925 ± 20 ± 46 ± 117) MeV , (22) 

in remarkable agreement (within the errors) with the experimental candidate X(3872) and with an 
estimate from relativistic quark model [13] 6 . The first error comes from sq, the second from the truncation 
of the OPE, and the third from the QCD inputs, such as m c and (qq). Despite this large dependence 
of our results on the value of the continuum threshold, the error induced by sq is comparable with the 
ones from other sources. The errors due to s$ could be reduced by using a more involved parametrization 
based on some effective Lagrangian or eventually, alternatives forms of the sum rules, such as those used 
in (48| for describing the hadronic r decay or some other sum rules [49| ■ 

Assuming that the mass of the first radial excitation is given by ^/sq, one can deduce a crude estimate 
of the splitting: 

X' - X w 225 MeV , (23) 

which is expected to be valid if the local quark-hadron duality is at work 7 . Within this assumption, one 
can notice that the mass-splitting is much smaller that the naive extrapolation from the ordinary meson 



6 However, due to the large error of our result, it can also be compatible with the X(3940), Y(3940) and Z(3930) if some 
of them are found to be a 1+ + state. We plan to study carefully the splitting of the different states of a given spin in a 
future work. 

7 If one uses similar assumption for the £>/(()+"'"), one can identify y^so with the radial excitation predicted in [30| | using 
some other approaches. 
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spectrum. Such a situation has been also encountered in the analysis of the pentaquark sum rule [43| 
and, in gen eral, in the analysis of correlators described by high-dimension operators such as hybrids and 
gluonia fl9j |. This result can indicate the existence of higher states near the lowest ground state mass, 
which can manifest as large continuum in the data analysis. 

One can also deduce to leading order in a s , from the individual lowest moments, the decay constant 
defined in Eq. ©: 

fx = (4.66 ± 0.16 ± 0.29 ± 0.68) x 10" 5 GeV , (24) 

which can be more affected by the radiative corrections than Mx- The first error comes from so, the 
second from the truncation of the OPE, and the third from the QCD inputs, such as m c and (qq). fx 
is useful for the estimate of its hadronic width using vertex sum rules. As X is an axial-vector meson, 
its decay constant can measure its weak transition into Iv via a VT-exchange,which might be difficult to 
measure experimentally, ft would be useful to have a measurement of this decay constant from some 
other methods, like e.g. lattice calculations. 



VI. SU(3) BREAKINGS AND MASS OF THE X 
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FIG. 6: The ratio Mx-> /Mx as a function of so, obtained from the LSR results up to different dimensions in the 
OPE. 

It is straightforward to extend the previous analysis to the case of the strange quark by using the QCD 
parameters given in Eq. (fTB|) . One can e.g. work with the ratios given in Eqs. ([Po]) and (|2"0)) . However, 
as the errors in the determination of Mx are relatively large, it will be difficult, to extract the SU(3) 
splitting from these individual ratio of moments. 

For extracting this relatively small mass-splitting, it is appropriate to use the double ratio of moments 

mil: 



M\ 



(25) 



for the LSR and for the FESR, which suppress different systematic errors (m c ,...) and the dependence on 
the sum rule parameters (sq,M 2 ). The results of the analysis from LSR are given in Fig. EJfrom which 
one can deduce, with a good accuracy: 



'dS. = 0.984 ± 0.002 ±0.007 



(26) 



where the first error comes from the QCD and sum rules parameters including the SU(3) breaking of the 
quark condensates. The second error from an estimate of the truncation of the OPE. This leads to the 
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mass splitting: 

M x > - M x ~ -(61 ± 30) MeV . (27) 

Similar methods used in [lj|[5(| have predicted successfully the values of Mb J Mb and Mb, /Mb, which 
is not quite surprising, as in the double ratios, all irrelevant sum rules systematics cancel out. 

Using FESR and taking the SU(3) breaking correction for the continuum threshold, which is the most 
important effect in the FESR analysis, we confirm the previous LSR result. It is interesting to notice that 
we predict a X s mass slightly lighter than the X, which is quite unusual. This is due to the fact that, 
in the sum rules expression of M Xs , the linear quark mass term tends to decrease the X s -mass, which 
is partly compensated by the effect of the quark condensates. Such a small and negative mass-splitting 
is rather striking and needs to be checked using alternative methods. Note, however, that a partial 
restoration of SU(3) symmetry is already observed in the neighbourhood of heavy quarks, illustrated 
by the almost equal hyperfine splittings D* s — D s and D* — D. In potential models, the mass spliting 
Mx° — Mx is certainly larger than the value in Eq. (j2"7|) , but smaller than 2(m s — m q ) as the increase of 
the constituent mass from m q to m s is partially cancelled by the deeper binding of the strange quarks. 
The existence of the X s , which can be experimentally checked, can serve for a further test of the four- 
quark model for the X. The (almost) degenerate value of the X and of the X s masses may suggest that 
the physically observed X state can result from a mixing between the ccqq and cess bare states, which 
may be dominated by its ccqq component. However, we expect that a careful and perfect analysis of the 
cess sector should feel the X in the spectrum, though with a small coupling. One should also notice that 
these ccqq and cess components can be comparable if the X is a SU (3) singlet state. 

Using the ratio of the s- over the q-quark sum rules , one can predict also the ratio of decay constants: 

^ ~ 1.025 ±0.010 (28) 
Jx 

where, in the individual sum rules, the m s corrections act positively implying that this ratio is larger 
than 1. We also expect the reliability of a such result advocating the previous arguments for the ratio of 
mass. Similar sum rule leads to Jb s /Jb = 1.16 ±0.03 (5l| . which has been confirmed later on by different 
lattice calculations. 

However, despite the different successful predictions of the ratio of moments for the i3-meson param- 
eters, we expect that the method will be less predictive for the four-quark state. This can be signaled 
by the large error in the previous prediction of the mass-difference. The inclusion of radiative or some 
other higher dimension condensates corrections or some other effects not accounted for in this paper will 
be useful for confirming or disproving the previous results. 



VII. TEST OF THE ISOSPIN VIOLATION 



We attempt to use of the sum rule, for a rough estimate of the small mass difference M(Xh) — M(Xi) 
defined in Eq. (TJJ. Using Eq. ([151) , we § et: 

M 2 (X h ) - M 2 (Xt) = 4m - So , (29) 

where: 

pi(s) = cos 2 9 p u (s) + sin 2 9 p d (s) and p h {s) = sin 2 9 p u (s) + cos 2 9 p d (s). (30) 

Here, p u {s) and p<i{s) are simply the spectral density p(s) defined before with the flavor of the light quark 
chosen as u and d, respectively. 

Clearly the only terms depending on the light quark flavor will contribute to the numerator of Eq. (129|) . 
In fact the expression ph(s) — pi(s) can be written in terms of the isospin breaking quantities: (uu) — (dd) = 
-j(qq), m u - m dj m u (uu) - m d (dd) = -^(qq)m q + (qq)(m u -m d ), and (to) 2 - (dd) 2 = -2^(qq) 2 , where 
(qq) = ((uu) + (dd))/2 and 7 = ((0\dd - uu\0))(0\uu\0). 

The value of 7 has been estimated in a variety of approaches with results varying over almost one order 
of magnitude: — 1 x 10~ 2 < 7 < —2 x 10 -3 [54(. However, studies based on chiral perturbation theory 
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[52| . LSR and FESR [53[ analysis of the (pseudo)scalar channels , analysis of the neutron-proton 
mass difference [54| and of the heavy meson decay widths [55[ leads to the value: 

7~ -(1±0.5) x 1CT 2 (31) 

1/2 

which we shall consider in our analysis . The results for the mass difference using s — 4.2 GeV can 
be seen in Fig. where we have considered two values for the mixing angle: 9 — 0° , corresponding to 
maximal isospin violations, and 6 = 20° which was the value determined in [l[ . 




M 2 (GeV 2 ) 

FIG. 7: M 2 (X h ) - M 2 (Xi) as a function of the sum rule parameter (M 2 ) for 7 = -1 X 1(T 2 . The solid line is for 
6 = 0° and the dashed one is for 9 = 20°. 

On can notice from Fig. [7] that the sign of the mass difference is reversed when one includes the 
dimension-six condensates, while the effect of the (partial) dimension-8 contribution is relatively small, 
indicating that the OPE starts to behave quite well. However, one needs a more complete evaluation 
of the dimension-8 contributions for a more precise determination of the mass-spiltting. For a more 
conservative estimate, we consider the range of the absolute value of the mass-difference, which is not 
strongly affected by the truncation of the OPE. In this way, one obtains from Fig. [7] 

\M{X h ) ~ M{Xi)\ ~ (2.6 - 3.9)MeV , (32) 

which is smaller than the (8± 3) MeV value given in [l[, but larger than the decay width of the A (3872), 
which is less than 2.3 MeV. However, one can notice from Fig. [7]that the sum rule cannot fix with a good 
precision the sign of the mass splitting, though it is tempting to conclude that the sign is negative, in 
disagreement with the result of [l[. 



VIII. SUM RULE PREDICTIONS FOR X b AND X£ 

Using the same interpolating field of Eq. (TT|) with the charm quark replaced by the bottom one, 
the analysis done for A(3872) in the previous sections can be repeated for Af,, where A& stands for a 
(bbqq) tetraquark axial meson. Using consistently the perturbative AfS'-mass mb(mb) = 4.24 GeV, and 
working with the LSR, we find a good OPE convergence for M 2 > 5 GeV 2 . We also find that, for 
so < (10.2 GeV) 2 , the continuum contribution is always bigger than the pole contribution for all values 
of M 2 > 5 GeV 2 . 

In Fig.[8]we show the Af, meson mass obtained from Eq. (|15|) . in the relevant sum rules window, with the 
upper and lower validity limits indicated. Although we get a good OPE convergence for M 2 > 5 GeV 2 , 
we have now a more restricted lower limit given by Mx b < \fso- Therefore, the lower limit indicated in 
Fig. [5] is given by this condition. 

From Fig. [8] we see that the results are very stable as a function of M 2 in the allowed region. However, 
the LSR prediction increases with sq. Taking into account the variation of M 2 and choosing (a priori) 
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9.5- 
9.0- 
8.5- 



o.u- 1 — | 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 

M 2 (GeV 2 ) 

FIG. 8: The Xt meson mass as a function of the sum rule parameter (M 2 ) for different values of the continuum 
threshold: y/IE = 10.2 GeV (solid line), = 10.5 GeV (dashed line) and ^ = 10.8 GeV (dotted line). The 
arrows delimit the region allowed for the LSR sum rules. 



some range of sq, we arrive at the predictions: 

10.06 GeV < M Xb < 10.50 GeV , 
for 10.2 GeV < s 1/2 < 10.8 GeV and 5.0 < M 2 < 8.5 GeV 2 . 

11.0- 

FESR (up to dim 5) 
LSR for ^=5.6 GeV (up to dim 5) 

10.5- 



> 

I 

sP10.0- 



(33) 




9.5- 



i 1 I 1 I 1 I 1 I 1 I 1 I 
10.2 10.3 10.4 10.5 10.6 10.7 10.8 

s 1/2 (GeV) 

FIG. 9: The FESR results in Eq. J20J for mx b as a function of s 1/2 for n = 0. The LSR result with M 2 = 5.6 GeV 2 
has been also inserted for easy comparison. 

The FESR analysis can also be repeated in the case of the 6-quarks for improving the LSR results. The 
results are shown in Fig. [5] As in the case of X(3872) the curves for n = and n = 1 are quite similar, 
and again there is no stability in sq. We also show in the same figure the LSR results for M 2 = 5.6 GeV 2 , 
as a function of sq 1 / 2 . A common solution is obtained for: 

sq 1 / 2 = (10400 ± 20) MeV , (34) 

to which corresponds the improved final prediction: 



M Xb = (10144 ±21 ± 104) MeV 



(35) 
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The second error comes again from the QCD inputs. The central value in Eq. (|3"3"|) is close to the mass 
of T(35'), and appreciably below the B*B threshold at about 10.6 GeV. For comparison, the molecular 
model predicts for A& a mass which is about 50 — 60 MeV below this threshold [2J , while a relativistic 
quark model without explicit (bb) clustering predicts a value of about 133 MeV below this threshold 
[47j |. It would also be interesting to have the (unquenched) lattice results for this state in order to test 
our QCD based results. A future discovery of this state, e.g. at LHCb, will certainly test the different 
theoretical models on this state and clarify, in the same time, the nature of the A (3872). 

One can also notice, by assuming that the mass of first radial excitation is about the value of ^/so'- 



X' - X 



X b « (225 - 250) MeV 



(36) 



For completeness, we predict the corresponding useful value of the decay constant to leading order in a s : 

fx, = (6.9 ~ 7.1) x 10- 6 GeV . (37) 

Our previous results will be useful inputs for studying more precisely the phenomenology of the Xf, 
outlined in 15611. 



1.04 
1.02 

1 



, 0.98 



' 0.96 ^ 



0.94 



0.92 



0.9 




10.2 



10.3 



10.4 



10.5 

' 2 (GeV) 



10.6 



10.7 



10.8 



FIG. 10: The ratio Mj< /Mx b as a function of so, obtained from the LSR results up to different dimensions in the 
OPE. The continuous curve corresponds to the D — 5 condensate contributions. The dotted and dashed curves 
correspond respectively to the vacuum saturation (K — 1) and violation of the vacuum saturation (K — 2) by a 
factor 2 of the dimension six- and (partial) eight-contributions. 

We extend the analysis to the X£(bbss). We show in Fig. [TO] the LSR prediction for the mass ratio, 
from which we deduce, by truncating the OPE at D = 5: 



M 



±0.002 ±0.018 



(38) 



where the first error comes from sq and the QCD parameters, while the second one from the truncation 
of the OPE. This leads to: 



M x ° - M Xb = -(123 ± 182) MeV , 



and 



fx b 



1.12 ±0.03 



(39) 



where the error due the truncation of the OPE is larger than in the case of the c-quark, which is mainly 
due to the terms of the form to^/M 2 in the OPE. We expect that the X^-family will show up at LHCb 
in the near future, which will serve as a test of our previous predictions. 



IX. CONCLUSIONS 



We have presented a QSSR analysis of the two-point functions of the X (3872) meson considered as a 
four quark state. We find that the sum rules result in Eq. (j2"2"|) is compatible with experimental data. An 
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improvement of this result needs an accurate determination of running mass m c of the MS'-scheme and 
the inclusion of radiative corrections. 

We have extended the analysis for studying the mass splitting between the X s and A due to SU(3) 
breaking. Our result in Eq. (|27| indicates an unusual ordering which deserves further independent checks 
from other QCD-based approaches, especially lattice calculations. However, our small mass-splitting 
suggests that perhaps the observed X has a cess component, though with a small coupling, with a size 
which depends on the SU(3) assignment of the X. 

Allowing possible isospin violations, we have also studied the mass splitting between the states Xi and 
Xh- Using the common values from different approaches of the leading SU(2) breaking parameter 7 of 
the light quark condensates defined in Eq. (l3"Tj) . we obtain the splitting in Eq. (f3"2"|) which is smaller 
than the value 8 MeV predicted in [l|. 

There are limits 58[ on the production of charged partners of the A (3872), but based on the weak 
decay of B mesons. It cannot be excluded that B decay favours neutral heavy mesons, if it proceeds first 
via an excited (cc) state recoiling against a cluster of light quarks and antiquarks. Hence the search for 
charged partners should be extended to other production mechanisms. 

Extending our analysis to the 6-quark meson, we found the values of the Af, and X£ masses in Eqs. 
([53)) and (j3U|) , which are appreciably below the B*B threshold at about 10.6 GeV. This is a common 
feature of all quark models with proper account for the correlation between the heavy quark and the 
heavy antiquark that the Af, is more deeply bound with respect to BB* than X c with respect to DD* , 
for the same reason why the (bb) family has more narrow states than (cc). In contrast, the molecular 
model, in which Af, is a meson-meson system bound by nuclear forces predicts this state rather close 
below the BB* threshold. Our analysis also indicates that the mass-splitting between the ground state 
and the first radial excitation is about (225 ~ 250) MeV, which is much smaller than the one expected 
from ordinary mesons, and which are (within the errors) flavour independent. 

We present in Eqs. (|24p . and lj3"T|) predictions of the decay constants of the X and Xb, and in Eqs. 
(I28p and (l3^|) the ratio of the strange over the non strange decay constants. These are useful quantities 
for studying the leptonic and hadronic decay widths of such mesons, and which can be checked from 
(unquenched) lattice calculations or from some other models. 

A future discovery at B factories or LHCb of the different states which we have predicted as a conse- 
quence of the 1 ++ four-quark nature assumption of the A" (3872) will certainly test the different theoretical 
models proposed for this state and clarify, in the same time, the nature of the A (3872). 

Different choices of the four-quark operators have been systematically presented for the ++ light 
mesons in [3l|, which should mix under renormalizations |32j from which one can deduce a "physical" 
renormalization group invariant current (RGI) which can describe the observed state. Though some 
combinations can provide a faster convergence of the OPE, we do not expect that the choice of the 
operators will affect much our results, where, in our analysis, the OPE has a good convergence while 
the renormalization mixing is a higher order effect in a s . Another choice of operator not included in 
the previous analysis is, for instance, given in ref . 0] , where it was shown that a simple chromomagnetic 
model suggests that the color octet-octet in the [cc]s=i [qq]s=i basis is the most natural candidate for 
describing the A (3872). However, this choice would correspond to an operator of higher dimension than 
the one analyzed in Eq. (fTJ), which would therefore induce relatively small corrections to the present 
analysis. Though done with a particular choice of current [l|, we expect that the results given in this 
paper will reproduce (within the errors of the approach) the general features of the four-quark model for 
the A (3872). We plan to come back to these issues in a future publication. 

Once the mass of the A (3872) is understood, it remains to explain why it is so narrow. There are 
presumably many multiquark states, but most of them are very broad and cannot be singled out from 
the continuum. In a recent study [57|, based on the same interpolating field as the one used here, it was 
shown that, in order to explain the small width of the A(3872), one has to choose a particular set of 
diagrams contributing to its decay. However, it will be desirable if this investigation can be checked from 
alternative approaches, such as lattice calculations. If confirmed, this method can be straightforwardly 
repeated to a variety of currents for understanding the width and the internal structure of the A (3872). 

If the A (3872) is a four-quark state, as our analysis suggests in answer to the question raised in the 
title, a four-quark structure probably holds for the states seen near 3940 MeV and 4260 MeV, on which 
more experimental information is still needed. This is our intention to extend the present analysis to 
other J PC configurations which are likely to host multiquark resonances. 



15 



Acknowledgements 

This work has been partly supported by the CNRS-IN2P3 within the (Non) Exotic Hadrons Working 
Group Program, and by FAPESP and CNPq-Brazil. S.N. wishes to thank the CERN-Theory Division 
where this work has been completed. 



[1] L. Maiani, F. Piccinini, A.D. Polosa, V. Riquer, Phys. Rev. D71, 014028 (2005). 

[2] E. S. Swanson, Phys. Rept. 429, 243 (2006); J. M. Richard, talk given at QCD 05 (Montpellier 4-8th July 

2005), hep-ph/0601043 
[3] S.-K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 91, 262001 (2003). 

[4] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 93, 162002 (2004); D. Acosta et al. [CDF II 
Collaboration], Phys. Rev. Lett. 93, 072001 (2004); B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 
71, 071103 (2005). 

[5] S.-L. Choi [Belle Collaboration], hep-ex/0405014 

[6] K. Abe et al. [Belle Collaboration] , |hep-ex/0505037[ |hep-ex/0505038| 

[7] G. Gokhroo and G. Majumder et al. [Belle Collaboration], arXiv:hep-ex/0606055 

[8] H. H0gaasen, J.-M. Richard, P. Sorba, Phys. Rev. D73, 054013 (2006); F. Buccella, H. H0gaasen, J.-M. 
Richard, P. Sorba, [hep^ph/0608001] Y. Cui et al. hep-ph/0607226 

[9] D. Bugg, Phys. Lett. B598, 8 (2004). 
[10] B.-A. Li, Phys. Lett. B605, 306 (2005). 
[11] K. K. Seth, P hys. Lett. B 612, 1 (2005). 
[12] N. Tornqvist, |hep-ph/0308277| 

[13] F.E. Close and P. Page, Phys. Lett. B578, 119 (2004). 

[14] C.Y. Wong, Phys. Rev. C69, 055202 (2004). 

[15] S. Pakvasa and M. Suzuki, Phys. Lett. B579, 67 (2004). 

[16] E. S. Swanson, Phys. Lett. B 588, 189 (2004); Phys. Lett. B 598, 197 (2004); [hep^ph/0601110 
[17] M.A. Shifman, A.I. and Vainshtein and V.I. Zakharov, Nucl. Phys. B147, 385 (1979). 
[18] L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127, 1 (1985). 

[19] For a review and references to original works, see e.g., S. Narison, QCD as a theory of hadrons, Cambridge 
Monogr. Part. Phys. Nucl. Phys. Cosmol. 17, 1-778 (2002) hep-h/0205006 ; QCD spectral sum rules , World 
Set. Lect. Notes Phys. 26, 1-527 (1989); Acta Phys. Pol. B26 (1995) 687; Riv. Nuov. Cim. 10N2 (1987) 1; 
Phys. Rept. 84,263 (1982). 

[20] R.A. Bertlmann, C.A. Dominguez, G. Launer and E. de Rafael, Nucl. Phys. B250, 61 (1985); R.A. Bertl- 

mann, G. Launer and E. de Rafael, Z. Phys. C 39, 231 (1988). 
[21] K. Chetyrkin, N.V. Krasnikov and A.N. Tavkhelidze, Phys. Lett. B76,83 (1978); N.V. Krasnikov, A. A. 

Pivovarov and A.N. Tavkhelidze, Z. Phys. C19,301 (1983). 
[22] J. Latorre and P. Pascual, J. Phys. Gil, L231 (1985). 
[23] S. Narison, Phys. Lett. B175, 88 (1986). 
[24] T.V. Brito et al, Phys. Lett. B608, 69 (2005). 
[25] H-J. Lee and N.I. Kochelev, hep-ph/0608188 
[26] M.E. Bracco et al., Phys. Lett. B624, 217 (2005). 

[27] H.-Y. Cheng, W.S. Hou, Phys.Lett. B566 (2003) 193; A. Datta, P. J. O'Donnell, Phys.Lett. B572 (2003) 164. 

P. Colangelo, F. De Fazio, R. Ferrandes, Mod. Phys. Lett. A19 (2004) 2083 ; A. Hayashigaki, K. Terasaki, 

hep-ph/0411285; D. Becirevic, S. Fajfer and S. Prelo vsek, |"hep-ph/0406296| Z.-G. Wang, Shao-Long Wan, 

hep-ph/0602080; J.P Pfannmoeller , hep-ph/0608213 
[28] S. Narison , |hep-ph/0512256| S. Narison, Phys. Lett. B605, 319 (2005); S. Narison, Phys. Lett. B210, 238 

(1988). 

[29] W.A. Bardeen et al., Phys. Rev. D68 (2003) 054024; S. Godfrey, Phys. Lett. B568 (2003) 254; M. Harada 

et al, Phys. Rev. D70 (2004) 074002; Y.B. Dai et al. Phys. Rev. D68 (2003) 114011. 
[30] O. Lakhina and E.S. Swanson, hep-ph/06080 lT| 
[31] H-X. Chen, A. Hosaka and S-L. Zhu, |hepq3h/ 0604049 

[32] S. Narison and R. Tarrach, Phys. Lett. B125, 217 (1983); M. Jamin and M. Kremer, Nucl. Phys. B277, 349 
(1986). 

[33] B. L. Ioffe, Nucl. Phys. B188, 317 (1981); B191, 591(E) (1981). 

[34] For a review, see e.g. H.G Dosch, Conference on Non-Perturbative Methods, Montpellier 1985, ed. World 
Scientific Co and references therein; H.G. Dosch, M. Jamin and S. Narison, Phys. Lett. B220 (1989) 251. 

[35] K. Chetyrkin, S. Narison and V.I. Zakharov , Nucl. Phys. B5 50 (1999) 353. 

[36] For recent reviews, see e.g.: V.A. Zakharov, hep-ph/0509114| S. Narison, [hep-ph/0508259| 
[37] S. Narison and V.I. Zakharov, Phys. Lett. B522 (2001) 266; K. Chetyrkin, S. Narison, Phys. Lett. B485 
(2000) 145; S. Narison, Phys. Lett. B300 (1993) 293; S. Narison B361 (1995) 121; S. Narison, Phys.Lett. 



16 



B626 (2005) 101. 

Bagan et al., Nucl. Phys. B254 (1985) 55; D.J. Broadhurst and S. Generalis, Phys. lett. B139 (1984) 85. 
S. Naris on, Nucl. Phys . Proc. Suppl. 86 (2000) 242 ( |hep-ph/99 114541 ); S. Narison, hep-ph/0202200 S. 
Narison, hep-ph/0510108; S. Narison, Phys. lett. B341 (1994) 73; H.G. Dosch and S. Narison, Phys. lett. 
B417 (1998) 173; S. Narison, Phys. lett. B216 (1989) 191. 
S. Narison, Phys. Lett. B466, 345 (1999). 

S. Narison, Phys. Lett. B361, 121 (1995); S. Narison, Phys. Lett. B387, 162 (1996). S. Narison, Phys.Lett. 
B624 (2005) 223 

A. G. Oganesian, hep-ph/0510327 

R. Matheus and S. Narison, Nucl. Phys (Proc. Suppl.) 152, 263 (2006) ( |hep-ph/0412062| ). 
G. Launer, S. Narison, R. Tarrach, Z. Phys. C26, 433 (1984). 

J. Bordes, V. Gimenez and J. A. Pcnarrocha, Phys. Lett. B201, 365 (1988). 
M. Davier, A. Hocker, Z. Zhang, hep-ph/0507078 and references therein. 

D. Ebert, R.N. Faustov and V.O. Galkin, Phys. Lett. B634, 214 (2006). 

E. Braaten, S. Narison and A. Pich, Nucl. Phys. B373 581 (1992); F. Le Diberder and A. Pich, Phys. Lett. 
B286 147 (1992). 

K. Maltman, Phys. Lett. B 440 367 (1998). 

S. Narison, Phys. Lett. B210, 238 (1988). 

S. Narison, Phys. Lett. B320, 369 (1994). 

J. Gasser, H. Leutwyler, Nucl. Phys. B250, 465 (1985). 

C. Dominguez and E. de Rafael, Ann. Phys. (N.Y.) 174 372 (1987). 

X. Jin, M. Nielsen, J. Pasupathy, Phys. Rev. D51, 3688 (1995); H. Forkel, M. Nielsen, Phys. Rev. D55, 1471 
(1997), and references therein. 

M. Nielsen Phys. Lett. B 634, 35 (2006) ( [hep-ph/05 102771 ). 
W.-S. Hou, hep-ph/0606016] 

F. S. Navarra, M. Nielsen, Phys. Lett. B639, 272 (2006) ( |hep-ph/0605038[). 

B. Aubert et al. [BaBar Collaboration], Phys. Rev. D71 (2005) 031501 flhep-ex/041205l} . 



